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Abstract 

Recent experimental and theoretical studies suggest that the quarkonia suppression in a thermal 
QCD medium created at heavy ion collisions is a complex interplay of various physical processes. 
In this article we put together most of these processes in a unified way to calculate the charmonium 
survival probability (nuclear modification factor) at energies available at relativistic heavy ion 
collider (RHIC) and large hadron collider (LHC) experiments. We have included shadowing as the 
dominant cold nuclear matter (CNM) effect. Further, gluo-dissociation and collision damping have 
been included which provide width to the spectral function of charmonia in a thermal medium 
and cause the dissociation of charmonium along with usual colour screening. We include the 
colour screening using our recently proposed modified Chu and Matsui model. Furthermore, we 
incorporate the recombination of uncorrelated charm and anti-charm quark for the regeneration 
of charmonium over the entire temporal evolution of QGP medium. Finally we do the feed-down 
correction from the excited states to calculate the survival probability of charmonium. We find 
that our unified model suitably describes the experimental nuclear modification data of J/f} at 
RHIC and LHC simultaneously. 
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I. INTRODUCTION 


The physical picture of quarkonium dissociation in a thermal medium has undergone 
theoretical and experimental refinements over the last decade |l|. Heavy quarkonia (J/V’, T 
etc.) suppression is considered as the most classical observable of QGP formation in heavy 
ion collision experiments. This is because the heavy mass scale (m = 3.1 GeV for J/ifj and 
m = 9.2 GeV for T) makes these system possible for analytical treatment theoretically. 
On the other side, decay of heavy quarkonia via dileptonic channel lead to relatively 
clean signal which can be precisely measured experimentally. Till the mid-2000s, Debye 
screening was thought to be the only possible mechanism for the anomalous suppression 
of charmonium (J/ip) and bottomonium (T) in QGP medium. However, experimental 
results involve some puzzling features which defy explanations based on color screening 
alone js-?]. The first such experimental result is the less suppression at mid-rapidity 
than forward rapidity observed at RHIG and also at LHG jb, 7| which is in contradiction 
to the color screening scenario (as color screening predicts larger suppression at higher 
density region of plasma which is actually the mid-rapidity). Second such experimental 
result is the same amount of charmonium suppression at SPS and RHIG energies for the 
same number of participants jd]. Although, the available energy spans over two order 
of magnitude in moving from SPS to LHG, the amount of charmonium suppression is 
found to be similar. Regeneration of charmonia in QGP through the recombination of 
c and c quarks is believed as the main reason for this experimental observation. Third 
experimental observation is the suppression pattern in forward and backward rapidity in 
d — Au collision at RHIG. A suppression is observed at forward rapidity (in the d-going 
direction) and an enhancement at backward rapidity (in the Au-going direction) js]. This 
result suggests the importance of charmonia break-up effects in nuclear matter at hnal 
stages of collision apart from usual cold-nuclear-matter effects in the initial stage. All these 
experimental observations suggest that the charmonium suppression in QGD plasma is 
not the result of a single mechanism, but is a complex interplay of various physical processes. 


On the theoretical side, the development in thermal field theory shows that the static 
potential between two heavy quarks placed in a QGD medium consists of two parts SQ. 


It is the hrst part which represents the standard time-independent Debye-screened potential 












(Earlier it was thought that the Debye-screened potential is only the dominant part in 
heavy-quark potential based on which one can understand the dissociation of quarkonia 
in QGP). The second part of potential, other than the standard Debye-screened part, is 
imaginary and in the limit of f oo, represents the thermal decay width induced by 
Landau damping of the low frequency gau ge helds that mediate interaction between two 
heavy quarks Qj. Later Brambilla et ah Hj] have shown that the apart from thermal width 
originates from the imaginary part of the gluon self energy, singlet to octet transition of 
heavy meson resonance due to the gluonic interaction also contribute to the decay width 
of quarkonia states. In a QGP, gluons can collide with a color-singlet heavy quarkonium 
leading to its dissociation [l2| . Dissociation by the absorption of a single gluon is allowed as 
the color-octet hnal state of a free quark and anti-quark can propagate in the colored QGP 
medium, in contrast to the color-less hadronic medium. One of the earlier treatment of the 
dissociation of heavy quarkonium by the absorption of a El gluon (where El is the lowest 
or the spin-orbital wave function of gluons) was carried out by Peskin and 


Bhanot 
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14| . Therefore, the dissociation of quarkonia does not happen only due to the 


Debye screening but it can occur by the gluonic dissociation and the collisional damping 
(Landau damping) as well. 


The production of charmonium in nucleon-nucleon collisions is also an involved pro¬ 


cess 


15| . The charmonia production process in elementary hadronic collision e.g., p + p 


collision, begins with the formation of a cc pair; this pair can then either lead to open 
charm production or subsequently bind to form a charmonium state (about 10% for 
all charmonia) 1^. The dominant high energy production mechanism for charmoina is 
gluon-gluon fusion. Based on the scales involved, the production process of charmonia 
is believed to be factorisable into two parts: a charm and anti-charm quark produced 
through nucleon-nucleon collision is a perturbative QGD process |^. However, the 
formation and evolution of this pair into a meson is governed by non-perturbative QGD. 
Hence, heavy-quarkonia provide a unique laboratory which can explore the interplay of 
perturbative and non-perturbative QGD effects. A variety of theoretical approaches has 


been proposed in t 
nucleon collisions 
approaches 
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le literature to calculate the heavy quarkonium production in nucleon- 
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20| are the two theoretical methods based on QCD which are being used 





























in most of the quarkonium production and suppression models. However, the heavy 
quarkonium production mechanism is still a topic of intense debate. 


The heavy quarkonia production processes, in the case of nucleus-nucleus collision, are 


signihcantly affected by the nuclear environment 


23|. These effects are known as cold 


nuclear matter (CNM) effects. In most of the literature, three CNM effects are considered. 
The first and dominant CNM effect in the case of quarkonium production is shadowing. 
Change in the parton distribution function in the nucleus which control the initial parton 
behaviour and strongly depends on the collisional kinematics ( in the small x region nuclear 
parton distribution function is clearly suppressed compared to that of nucleon) is known 
as shadowing [2^. The shadowing cause the production cross-section to become less in 
A — A case to that of pure N-N collision. Second CNM contribution is known as Cronin 


effect 


25 


26[ |. It describes the initial gluon multi-scattering with the neighbouring nucleons 


presented in the nucleus prior to hard scattering and the quarkonia formation. This 
results in the broadening of transverse momentum distribution of produced charmonia. 
Nuclear absorption 0 is another CNM contribution to the charmonia production. The 
interaction between charmonium and the primary nucleons leads to the normal suppression 
of charmonia. Nuclear absorption is the dominant CNM effect at lower energies. However, 
the cross-section for nuclear absorption decreases with the increase in energy 281. 


Recently, it has been proposed that recombination o 


in QGP can also regenerate the charmonia states 29 


initially uncorrelated c and c quarks 


351] • The calculation of regeneration 


of charmonium is based on statistical hadronization model 


29 


30| and kinetic model in 


which the J/V’ production is described via c 


whole temporal evolution of the QGP 


31 33 


ynamical melting and regeneration over the 


35|. Some transport calculations were also 


performed to calculate the number of regenerated J/V's [Sg, l37|. At lower energies, this 
contribution is very low almost negligible because of the fewer number of initially produced 
charm quarks. However, at higher energies of RHIG and LHG, the regeneration factor 
become important. Thus J/V’ whose suppression is actually proposed as a signal to confirm 
the existence of quark gluon plasma earlier can also turn out to provide extremely useful 
probe for QGD medium created at heavy ion collision experiments. 








































In this article we present a unified model which includes most of the above discussed 
dissociation as well as production (recombination) processes to finally calculate the survival 


probability of J/'0 in QGP medium. We have constructed this model 


Dased on the kinetic 


approach whose original ingredients was given by Thews et ah 33l-l35| . In this approach, 


there are two terms written on the basis of Boltzmann kinetic equation. First term, which 
we call as dissociation term, includes the dissociation processes like gluonic dissociation and 
collisional damping. The second term (formation term) provides the (re)generation oi J/'ip 
due to the recombination of charm-anticharm quark. These two terms compete over the 
entire temporal evolution of the QGP and at freeze-out temperature we get the multiplicity 
of finally survived J/ijjs. To include the gluo nic dissociation we took help of a model which 


was developed by Nendzig and Wolschin 


and later used by two of the authors |39(] . The 


thermal width due to collisional damping is calculated based on the thermal field theory as 
discussed by M. Paine et. ah jo]. We have also included the shadowing effect to incorporate 
the GNM effect properly into the production process. We consider color screening as a 
dissociation process of charmonium acting till the formation of charmonium bound states 
followed by gluonic dissociation along with the collisional damping. It means that the 
color-screening is active at the initial times of medium evolution when the temperature is 
high enough to melt-down the charmonia states. Later the dissociation probability by color 
screening diminishes rapidly and becomes zero at lower temperatures (time larger than the 
quarkonia formation time). To include the dissociation oi J/pj due to Debye screening (color 
screening) in the QGD plasma, we have used a new model constructed by two of the authors 
based on the color screening in the QGP 411. In this color screening model we have used 
the quasipaticle (QPM) equation of state (EOS) to describe the basic partonic properties of 
QGP phase. To define the dynamics of the system created in the heavy ion collisions, we 
have used the (l-l-l)-dimensional viscous hydrodynamics. We have included only the shear 
viscosity and neglect the bulk viscosity. We have also suitably incorporated the overall 
feed-down correction from the higher charmonium states (xc and ip') to the charmonium 
(J/V^). 
















II. MODEL FORMALISM 


The abundance of charm quark, anti-quark and their bound state i.e, charmonia states 
{J/'ip, Xc, 4’' etc.) is governed by a simple master equation involving two reactions: the 
formation reaction and the dissociation reaction. Thus the time evolution of the number of 


bound charmonium state in the deconfined region can be written as 33l |: 

dNj/^ 


dr 


= N, [V{t)]-^ - Tn,niNj/^. 


( 1 ) 


In the above equation, the first term in the right hand side represents the formation term 
by recombination of uncorrelated charm quark and anti-quark. The second term in the 
right hand side is the dissociation term of charmonium. Tand T p^ni are the dissociation 
width and recombination reactivity corresponding to the dissociation and regeneration of 
charmonia, respectively. It is important here to mention that the unit of Tis in GeV 
or fm~^. However, the recombination reactivity Tp^ni has its units in frr4 or GeV“^. It 
only change in the unit of fm~^ or GeV when it is multiplied by the inverse of the system 
volume [H(r)]“^. Nc, Nc and Nj/^ are the numbers of produced charm, anti-charm and 
J/V', respectively. At the initial time, we have taken Nc=Nc=Ncc- When the total number 
of regenerated charmonia is less than the initial number of N^c, one can obtain the analytical 
solution of Eq. ([1]) as follows [3^ : 


Nj/^{Tf,b) = e{Tf) 


fTf 


A^j/p(t‘o,&) + / rir,„i[l/(r)e(r)] ^dr 


' TO 


( 2 ) 


where Tq is the initial thermalization time of the QGP and tj is the life time of the QGP. 
Nj/^{To,b) is the initial multiplicity and Nj/^{Tf,b) is the hnally survived number of J/4 
meson. The variables used in the Eq. ([2]), e(r/) and e(r) are the suppression factors which 
can be obtained using the following expressions : 

rs 


and 


:(r/) = exp 


e[T) = exp 


D,nl 


dr 


' TO 


- / P 


D,nl 


dr 


(3) 


(4) 


We have dehned Tp^ni as the net sum of collisional damping reaction rate (J'damp,ni) and 
gluonic dissociation reaction rate (Pgd,nz) of charmonia in QGP, given as; 


^ D,nl — r 


gd^nl 


+ P 


damp,nl • 


(5) 











The initial time (tq) is taken as the formation time required for the charmonia formation 
and where the dissociation due to color screening becomes zero, tq is taken as 0.89, 2.0 and 
1.5 for J/'0, Xc and if}', respectively 421. 


A. Gluonic Dissociation and Collisional Damping 


The dissociation of cc bound state due to gluonic dissociation along with 
damping in QGP medium was formulated by Wolschin et el. Q H- Here 
discussed about these dissociation mechanisms. 


collisional 
we briefly 


1. Collisional Damping 


To determine the collisional dissociation and gluonic dissociation rate, we take the help 
of effective potential models. In our work, we have used the singlet potential for c — c bound 
state in the QGP as follows [38|; 

p-rrior' 
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( 6 ) 


where the hrst and second term in the right hand side is the string term and the columbic 
term, respectively. The third term in the right hand side is the imaginary part of the 
heavy-quark potential. In the above Eq. (6) : 

• cr is the string tension constant between cc bound state, given as cr = 0.192 GeV^. 

• is Debye mass, rriD = T H coupling constant at hard 


scale, as it should be aj = Q;s(27rT) < 0.50 


. For charmonia we found aT ~ 0.25 


. Here we have considered, Nc = 3&:Nf = 3 and evolution of temperature, T as the 
function of time, r and impact parameter, b 391. 


Qfe// is effective coupling constant, depending on the strong coupling constant at soft 
scale a® = as{mcas/2) = 0.4725, given as Oe// = 

















The decay rate, Tdamp,nh accounts for collisional damping by the QGP partons. The imag¬ 
inary part of potential causes the collisional damping (also termed as Landau damping in 
the literatures). Therefore, the decay rate can be obtained, in a hrst order perturbation, by 
folding of imaginary part of the potential with the radial wave function and is given by; 

r.„„„, = /b„,(r)M/m(r)|<,„,(OI<ir. ( 7 ) 

here, Qnii'r') is charmonia wave function. Corresponding to different values of n and I, we 
have obtained the wave functions for 1S'( lP{xc) and 2S{'ip') by solving the Schrodinger 
equation. 


2. Gluonic Dissociation 


In a QGP, gluons can collide with a color-singlet heavy quarkonium leading to its dis¬ 
sociation. The ultra-soft gluon makes the color-singlet object to a color-octet object which 
further dissociates at the time of freeze-out. The cross-section for this gluonic dissociation 
process can be given as follows 


ITP \ — 

^d,nl[Dgj — ^2 


TUr 
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nl 


Eg + Enl 


21 + 1 


( 8 ) 


where, is the probability density obtained using the singlet and octet wave functions as. 


rqV _ 


drr g*^i{r) hqv{r) 


(9) 


and 


rric = 1.3 GeV, is the mass of the c and c. 


~ 0.59 [38|, is coupling constant, scaled as = as(asm^/2). 


Enl is energy eigen values corresponding to charmonia wave function, gniir). 


• the radial wave function hqi>{r) has been obtained by solving the Schrodinger equation 
with the octet potential Vg = aeff/Sr and the value of q is determined from energy 
conservation, q = ^mdEg -|- Eni). 










FIG. 1: (Color online) Variation of total dissociation width with temperature along with its 
components i.e., gluonic dissociation width and width due to collisional damping. 


The Schrodinger equation has been solved by taking a 10^ point logarithmically spaced 


hnite spatial grid and solving the resulting matrix eigen value equations 


^|. For the 


octet modeling the potential is repulsive , which implies that the quark and antiquark can 
be far away from each other. To account for this, the hnite spatial grid is taken over a 
very large distance, namely 10^, as an approximation for inhnity. The octet wave function 
corresponding to large cc distance have negligible contribution to the gluonic dissociation 
cross section. 


To obtain the gluonic decay rate we have calculated the mean of gluonic dissocia¬ 

tion cross section by taking its thermal average over Bose-Einstein distribution function for 
gluons. Thus the rate of gluonic dissociation can be written as: 

P gd r dpgPlad,ni{Eg) 

gd,nl 27^2 _ 1 

where = 16 is the number of gluonic degrees of freedom. The Eq. ffTUD has been derived 
for idealized case where J/'0 is at rest in a thermal bath of gluons. 













FIG. 2: (Color online) Variation of gluo-dissociation cross section of J/ij; with respect to gluon 
energy, Eg at T = 0.170 GeV and T = 0.300 GeV. 


The total dissociation rate along with its components i.e., gluonic dissociation and colli- 
sional damping decay rates versus temperature are shown in Fig. 1. This figure depicts that 
the decay rate of charmonia due to gluonic dissociation remains insignificant till T = 0.250 
GeV and increases slowly beyond this temperature. Whereas, collisional damping decay 
rate contributes significantly throughout whole temperature range. Gluonic dissociation 
cross section for versus gluon energy, Eg at the temperatures, T = 0.170 GeV and 
T = 0.300 GeV are also depicted in Fig. 2. 


B. Regeneration via c and c quarks 


The recombination reactivity, Frequired in Eq. (2) is calculated by taking the thermal 
average of product of recombination cross section and initial relative velocity between c and 
c, < af^ni Vrei >p^ usiug modified Fermi-Dirac distribution function of charm quark at 


temperature T as follows 


35|; 


^ F,nl — 




Pc,m 


pI fc fc Vrel 




Pc,max CPc,max 


, r~dp,dp,pipihf, 

yc,rmn ^ yc^min 


( 11 ) 


where, Pc and Pc are momentum of charm and anti-charm quark, respectively. The /c,c 


is the modihed Fermi-Dirac distribution function of charm and anticharm quark given as, 
/c,c = Ac,c/-|-1), here = \/Pc,c + "^c,c is the energy of charm and anticharm quark 










in medium with mass, ruc^g = 1-3 GeV and Ac,g is their respective fugacity term 
have calculated relative velocity of cc pair in medium, given as: 


43j. We 


'^rel 


(pc - "i; 


( 12 ) 


pI pg + rnl{pl + Pg + ml) 

The recombination cross section Cf^ni has been obtained using detailed balance from the 
dissociation cross section ad,ni [33| as follows: 


48 , 


^f,nl / 

36 s{s 


4 ml) 


(13) 


Here, Mni = Xa i’'): is the mass of the charmonium states. The s = (pi + P 2 )^ is 

related to the center of mass energy of the cc pair with pi and p 2 as four momentum of c 
and c quarks, respectively. 


The Tp for J/ijj versus temperature T at both, RHIC and LHC energies are shown 
in Fig. 3. Its value at RHIC energy is found to be larger as compared to corresponding 
number at LHC energy. This trend may be due to the low initial momenta of charm and 
anti-charm quark which participate in the J/psi formation at RHIC compared to LHC. It is 
also obvious from the above hgure that the recombination reactivity. Tig at LHC increases 
with temperature, attains maximum value at around T = 0.225 GeV and vanishes at 
T = 0.300 GeV. Whereas, at RHIC energy it peaks at T = 0.300 GeV and remains hnite 
even at temperature T = 0.400 GeV. 


C. Inputs to the model 

In this section we provide the prescriptions to calculate various inputs which have been 
used in our model. Nj/^{to) and Ncc in Eq. (2) is the initially produced J/ip and cc pair 
in heavy ion collisions. We have calculated these quantities using Glauber model, for per 
event as follows; 

Nj/4To,b) = aJ/^TAA{b) (14) 


here, TAA{b) is nuclear overlap function, its impact parameter (6) dependent values have 
been taken from Ref. 4^. The is the cross section for production in p-|-p collision and 









FIG. 3: (Color online) Variation of recombination reactivity with respect to temperature at RHIC 
{y'sjsfN = 200 GeV) and at LHC {y/sjyN = 2.76 TeV). 


its values for RHIC at ^/s = 200 GeV and for LHC at ^/s = 2.76 TeV are given in Table I. 


Table 1. All cross sections are in mili-barn (mb) unit 


NN 

NN 

Xc 

-T 

LHC 

|— 1 

0.0072 45| 


o.3o</j: 

RHIC 

0.00139 

451 

i Orr^^ 

o-3o</i: 


Similarly, we have calculated from Glauber model; 


JVcc(l>) = ffS" Taa, 


( 16 ) 


where is the cross section for cc pair production in p+p collision. The has been 
calculated using pQCD approach for GRV HO hadronic structure function |3J], we have 
obtained = 3.546 mb for LHC at y/s = 2.76 TeV and = 0.346 mb for RHIC at 
y/s = 200 GeV. The quantity V (r) is the volume as a function of time r. It is based on the 
QPM EOS of QGP and the isentropic evolution of QGP ^ and given by, 


V(r, b) = vo{b) 


(16) 


where, R is the Reynolds number 40|] and vo{b) is the initial volume at time tq, given as, 
Vo{b) = 71 {rt — 6/2)^ro, here is the radius of hreball created. 































Here we use a cooling law for temperature which not only depends on proper time (r) but 
also varying with respect to number of participants {Npart)- The cooling law for temperature 
which connects proper time and Npart to the temperature of the system is as follows 39| : 

V3 1/3 

(17) 


rpt ( Npart{hin) \ ( Tqgp N 

^ ^ ^\Npart{Uno)) V r ) 


where Npart{bino) is the number of participant corresponding to the most central bin as used 
in our calculation and Npart{bin) is the number of participant corresponding to the bin at 
which we want to calculate the temperature, tqgp is the lifetime of QGP. 


D. Cold Nuclear Matter Effect 

We have already discussed about shadowing, absorption and cronin effect as the three 
main nuclear effects on the charmonium production. Nuclear absorption and Cronin effect 
are not included in our calculation. We incorporate shadowing as the only dominant CNM 
effect in the current work. 


1. Shadowing Effect 


We have used the EPS09 parametrization to obtain the shadowing for nuclei, with atomic 


mass number A, momentum fraction x, and scale /r, S'*(y4, x, fi) [4g, The spatial 


variation of shadowing can be given in terms of shadowing and the nucleon density z) 
as follows: 


Sl,{A, X, /i, r, z) = 1 + Np[S\A, x, /i) - 1] 


/ dz PA{r,z) 


f dz Pa 


where Np is determined by the following normalization condition: 


a 




39|, 


— J d^rdz pa{s) Sp{A, x, p,r, z) = x, p). 

The suppression factor due to shadowing is defined as; 

daAA/dy 


Ssh = RAA{b) = 


As mentioned in Ref. 


Taa d(Tppldy 


(18) 


(19) 


( 20 ) 


48|, the color evaporation model gives ct^a and as follows: 


<xaa = / dzi dz 2 d‘^r dxi dx 2 [fAA, xi, /i, r, zi) 


xfgiA X2, /i, b-r, Z2) agg^QQ{xi, X2, p)], 


(21) 












(22) 


^PP j dxi dX2 [fg{p, Xi, /i) fg(p, X2, p) 

^99—^ 2 ) /^)]- 


Here, xi and X 2 are the momentum fraction of the gluons in the two nuclei and they are 
related to the rapidity [39[. The superscripts i and j refer to the projectile and target nuclei, 
respectively. 

The function fg{A, x, p, r, zi) is determined from the gluon distribution function for 
proton fg{p, X, p) by using the following relations: 


• fg{A, xi, /i, r, zi) = pA{s)S\A,xi,p,r,z) fg{p, Xi, p). 

• fiiA, X 2 , p, b-r, Z 2 ) = pA{s)S^{A,X 2 ,p,b-r,z)fg{p, X 2 , p). 


The value of the gluon distribution function fg{p, x, p) in a proton (indicated by label p) 
has been estimated by using CTEQ 6 [49 1. 


E. Color Screening 


We treat color screening as independent suppression mechanism acting till formation 
of charmonium bound states followed by gluonic dissociation along with the collisional 


damping. Original color screening mechanism 


50| have been modihed by Mishra et ah, 51 |. 


In our present work we have used^asiparticle model (QPM) equation of state (EOS) based 


model in which pressure prohle 


40| | and cooling law of pressure are the main ingredients. 


The cooling law of pressure is given by 


. ^ . BCD 

p(r,r) = H+ — + — + ^ 


r<i 




(23) 


where A = -ci, B = C 2 C^, C = and D = here Ci, C 2 , C 3 are constants and have 


been calculated 


40 


4l| using different boundary conditions on energy density and pressure. 


Determining the pressure prohle at initial time r = and at screening time r = we get, 

( 24 ) 


BCD 

p{Ti, r) = A + — + — + — = p{Ti, 0) h{r) 













/ N . BCD 
p{Ts, r) = A + ^ + — + ^= pqgp, 


Q ■ ■ 

Ts Ts rt” 


(25) 


here Pqgp is the pressure of QGP inside screening region, required to dissociate J/'?/’, as 
determined by QPM EOS of QGP ^|. After combining cooling law and pressure profile 
and equating screening time to the dilated formation time, we determined the radius of 
screening region r*. 

Assuming that cc is formed inside screening region at a point whose position vector is 
f. It moves with transverse momentum pp making an azimuthal angle (j) (angle between 
the transverse momentum and position vector rj/^p). Then the condition for escape of cc 
without forming charmonium states is expressed as: 


cos (j) >Y; Y = 


(26) 


2rj/p,TFPT 

where, rj/^ is the position vector at which the charm, anti charm quark pair is formed, Tp is 
the proper formation time required for the formation of bound states of cc from correlated 
cc pair and m is the mass of charmonia (m = Mjjp;, for different resonance 

states of charmonium). 


1. Survival Probability 


In the color screening scenario, the survival probability of charmonia in QGP medium can 
be expressed as; 


2{a + 1) 

Sc{pT 1 ^part) ZPrP? / dr T (j^maxY) 


Jq 


.,2 'I “ 


1 - 


Rl 


where a = 0.5 


50 


m- 


The condition on azimuthal angle (j)max given by Eq.(4) is expressed as; 

TT if y < —1 


r = 


TV — cos ^|y| if 0 > y > — 1 
cos~^\Y\ if 0 < y < —1 


(27) 


0 if F > 1 














Then we have obtained pt integrated snrvival probability in the color screening scenario, 
given as [40|, ; 

Sc{N,art) = ( 28 ) 


rPTma. ^ 


for J/'0, Xc and -0' denoted by Sc^'^, S}'' and Sf , respectively. The range of px valnes nsed 
here is allowed by the corresponding experimental data. 


F. Net Survival Probability 


N)/^{TQ,h) = Nj/^{To,b) Ssh 

In Eq. ([2]), we have replaced Nj/^{to, b) with b) from Eq. 

Eq. ([2]), as follows; 

='(v) 


(29) 

and re-obtained the 


r-TJ- 


6) + Afi / Tf,,, |\/(T)£(T)]-‘<iT 


'TO 


(30) 


Now the snrvival probability dne to shadowing and glnonic dissociation along with collisional 
damping can be written as; 

NljTr.b) 

(31) 


oJ/'lp _ j/'tp^ J ' / 

^ Nj/i,{To,b) 


In this model we have assnmed that at initial stage of QGP the color screening would be 
dominated independent to the glnonic dissociation but at later stage it would partially affect 
the J/'0 formation and can be coupled with glnonic dissociation. So we have incorporated 
the color screening at the initial stage of QGP followed by glnonic dissociation along with 
the collisional damping including shadowing. We expressed the survival probability using 
Eq. fl28ll and Eq. (EH); 

gj/^ _ gjji> gm (32) 

In the same way we calculated the survival probability for and ip', written as, and 
S'j , respectively. 

It has been observed that only 60% oi J/ip come up by direct production whereas 30% is 
from the decay of Xc and 10% is form the decay of ip', so the net survival probability Sj/^ 
of a mixed system after incorporating feed-down is given as: 

0.60 Nj/^ + 0.30 Sf + 0.10 sf 


Sj/tjj — 


0.60 Nj/^ + 0.30 + 0.10 


(33) 










III. RESULTS AND DISCUSSIONS 


Before going to our results, we want to state that some short-hand notations have been 
used by us to show the different physical processes in the plots. We have used GD for gluonic 
dissociation, CD for collisional damping, CS for colour screening, S for shadowing as a CNM 
effect, FD for the feed-down correction and R for the regeneration via recombination of c 
and c quarks. 
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FIG. 4: (Color online) Survival probability (Sp) of J/ip in QGP medium versus Npart including glu¬ 
onic dissociation (GD) and collisional damping (CD) at RHIC energy with and without shadowing 
effect (S). 
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FIG. 5: (Color online) Survival probability versus Npart with gluonic dissociation (GD) and colli¬ 
sional damping (CD) along with color screening (CS) at RHIC with and without shadowing effect 
(S). 
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FIG. 6: (Color online) Same as Fig. 6 but recombination effect (R) is also included. Only recom¬ 
bination (R) is also shown here. 
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FIG. 7: (Color online) Same as Fig. 7 but feed-down (FD) due to higher resonances namely, Xc 
and is also included. 

Fig. 4 depicts our prediction on suppression (in terms of pt integrated survival probabil¬ 
ity, Sp) of JI'll) in the QGP medium formed at RHIC as a function of Npart at mid rapidity 
arising due to gluonic dissociation and collisional damping with and without shadowing 
as a CNM effect. The experimental data on suppression obtained from the PHENIX 
experiment at RHIC at center of mass energy, ^/snn = 200 CeV are also shown on 
the same graph for comparison with our predicted results. It is obvious from the above 
hgure that our result under-predicts the observed J/'ip suppression data without including 
shadowing. However, after including shadowing, it captures the experimental data on sup¬ 
pression reasonably well. In Fig. 5, we have plotted survival probability versus Npart after 
including color screening along with the gluonic dissociation and collisional damping with 



























and without shadowing effect. Our predicted result is almost similar to the shown in Fig. 4 
since color screening only affects on higher charmonium states at RHIC energy. In Fig. 6, 
we have plotted variation of recombination factor at RHIC energy versus Npart- This graph 
indicates that at RHIC energy, recombination for J/'^ increases very slowly with Npart and 
reaches to slightly greater than unity at the most central collisions. All the above mentioned 
suppression (GD+CD+CS+S) along with the recombination are also plotted on the same 
graph which obviously will be nearly identical to the variation shown in Fig. 5 since the 
recombination factor turns out to be small even with higher recombination reactivity at 
RHIC as compared to LHC. This trend of reactivity is due to the low initial momentum of 
c and c quarks (in contrary to that happens at LHC due to high initial momentum) which 
participate in the secondary charmonia formation at RHIC energy. Furthermore, the recom¬ 
bination is mainly governed by the number of cc pair produced initially and reactivity of the 
uncorrelated c and c pair. Due to less number of cc pair produced during the initial stage of 
collisions at RHIC energy and its quadratic dependence on Ncc, the recombination is found 
to be small even with signihcant recombination reactivity at RHIC energy. So far, we have 
not included feed-down due to the decay of higher resonances of charmonium namely, Xc and 
xfj' to J/'0- Therefore, survival probability of J/ip with contributions coming from all the 
suppression mechanisms (GD-|-CD-|-CS-|-S) mentioned above along with the recombination 
and feed-down via the decay of Xc and tjj' to J/'ifj is plotted with respect to Npart in Fig. 
7. After including the feed-down due to higher charmonium reasonances, J/ijj suppression 
increases at all centralities at RHIC energy but still within experimental error bars. Thus, 
our results agree reasonably well with the suppression data at mid rapidity obtained from 
the PHENIX experiment at RHIC at center of mass energy, y/s^N = 200 GeV. 

Fig. 8 shows J/'ip suppression at mid rapidity versus Npart due to the gluonic dissociation 
and collisional damping with and without shadowing as a GNM effect. The experimental 
data on J/'ll) suppression at mid rapidity obtained from the GMS experiment at LHG ^ 
center of mass energy, = 2.76 TeV are also depicted on the same graph for compari¬ 

son. Our result without shadowing shows good agreement with the experimental data while 
with shadowing effect, it over predicts the suppression. Fig. 9 includes color screening with 
the above suppression contributions namely, gluonic dissociation and collisional damping 
with and without shadowing effect. This figure implies that our predicted result captures 
data without inclusion of shadowing effect. While with shadowing effect, it shows over 
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FIG. 8: (Color online) Survival probability (Sp) of J/V' in QGP medium versus Npart including 
gluonic dissociation (GD) and collisional damping (CD) at LHC energy with and without shadowing 
effect (S). 
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FIG. 9: (Color online) Survival probability versus Npart with gluonic dissociation (GD) and colli¬ 
sional damping (CD) along with color screening (CS) at LHC with and without shadowing effect 
(S). 

suppression. Regeneration factor due to recombination of uncorrelated c and c pair at LHC 
energy and feed-down via decay of higher resonances are two important phenomenon which 
need to be incorporated to explain the data. Therefore, in the next two hgures, we have 
plotted above recombination effect with all the above mentioned suppression mechanisms 
(GD-I-CD-I-CS-I-S) with and without feed-down arising due to higher resonance states of 
the charmonium. Fig. 10 presents all the above suppression contributions and regeneration 
factor without including feed-down effect due to higher resonances. Only recombination is 
also presented on the same plot. This hgure indicates that recombination of uncorrelated cc 
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FIG. 10: (Color online) Same as Fig. 9 but recombination effect (R) is also included. Only 
recombination (R) is also shown here. 
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FIG. 11: (Color online) Same as Fig. 10 but feed-down (FD) due to higher resonances namely, Xc 
and ifj' is also included. 

pairs at LHC energy varies from around 1.1 at the most peripheral to around 1.3 at the most 
central collisions. Compratively larger recombination occurs here due to the sizable number 
of cc pairs produced at the LHC even with lower recombination reactivity. This hgure 
clearly indicates that our predicted results capture the trend of data spanning over almost 
the whole range of Npart- However, it depicts slight over suppression almost at all centralities. 


From Fig. 10, it is also evident that at LHC energy, recombination begins to decrease 
beyond a certain value of Npart in contrary to the expectation based on its quadratic de¬ 
pendence on Ncc pairs and on volume V{t). In fact, recombination has somewhat complex 
dependence since it does not only depend on the Ncc and V (r) but also on the recombina- 























tion reactivity Fof uncorrelated c and c which further depends on the temperature of the 
medium and momentum of c and c. Decrease of recombination factor at the highest Npart 
i.e., at the most central collision may be due to the peak value of F^? occurring at compra- 
tively lower temperature (than at RHIC energy) and signihcantly small value at the highest 
Npart- The temperature corresponding to that Npart at LHC is too large and at that much 
higher temperature, the value of F^? becomes very small (as can be seen from Fig. 3). That 
is why even with the highest value of N^c at the most central collisions at LHC energy, the 
small and the reducing trend of F^? with temperature gives the overall decreasing behavior 
to the recombination factor. One can observe from Fig. 6 that recombination factor at 
RHIC does not depict this kind of diminishing trend at the most central collisions. This is 
because the corresponding F p value remains signihcantly large at RHIC (see Fig. 3) even 
at the most central collisions (at high temperatures). Thus there is no reducing trend in the 
recombination factor at the most central collisions at RHIC unlike at LHC energy. Despite 
of the above facts, recombination factor at LHC is always higher than the corresponding 
value at RHIC energy for each centrality class (comparison between Figs. 6 and 10) which 
shows the dominance of dependence. 

Fig. 11 depicts our results for survival probability oi J/'ip with respect to centrality after 
including contributions from all the suppression mechanisms (GD+CD+CS+S) as a function 
of Npart along with the recombination and feed-down due to decay of higher charmonium 
resonance states to J/ip. Comparison of Figs. 10 and 11 shows that feed-down due to the 
decay of higher charmonium states io J/p) increases the suppression a little bit at LHC energy 
and slightly over predicts the suppression. However, it still depicts reasonable agreement 
with the data at LHC energy under the experimental uncertainties. Thus, our analysis 
shows that the current unihed model approach based on the combination of commonly 
employed suppression and recombination effects present reasonably good agreement with 
the experimental data at RHIC and LHC energy over the whole range of centrality. 


IV. CONCLUSIONS 


In conclusions, we have explained recent J/ip suppression data at mid rapidity obtained 
from RHIC and LHC experiments using a common formulation based on Debye color screen¬ 
ing, gluonic dissociation and collisional damping (popularly termed as Landau damping in 



the literatures) along with the recombination of uncorrelated c and c pairs in the later stage 
QGP formation. Shadowing of parton distribution as a CNM effect determined by Vogt 
approach and feed-down due to decay of higher resonances of charmonium have also been 
included in the current work. Recombination effect is incorporated via a transport approach 
using two self coupled transport equations. The QGP medium is assumed to be expanding 
under Bjorken’s scaling law at mid rapidity. Our current model based on the combination 
of color screening, gluonic dissociation along with collisional damping plus recombination 
effect explain J/'^ suppression data obtained from both the energies reasonably well without 
introducing any extra parameter while moving from one set of data at one energy to another 
set of data at other energy. The explanation of two different sets of suppression data at two 
different energies, differing by few orders of magnitude, employing a single set of mechanisms 
without including any new parameter, obviously add importance to our current approach. 
It is also clear from our results that J/ijj suppression in QGP medium is a result of many 
complex suppression mechanisms contrary to a single mechanism. Recombination is also 
important to include to explain charmonium suppression especially at LHG energy. It is 
worthwhile to note here that although there are few parameters in our current formulation 
yet not even a single parameter is varied freely in order to explain the suppression data. We 
have taken all parameter values based on earlier works. 
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